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Abstract. Previous studies of strong gravitational lensing by galaxy clusters have neglected the potential impact of the intra- 
cluster gas. Here, we compare simulations of strong cluster lensing including gas physics at increasing levels of complexity, 
i.e. with adiabatic, cooling, star-forming, feedback-receiving, and thermally conducting gas, along with different implementa- 
tions of the artificial viscosity in the SPH simulations. Each cluster was simulated starting from the same initial conditions so 
as to allow direct comparison of the simulated clusters. 

We compare the clusters' shapes, dynamics, and density profiles and then study their strong-lensing cross sections computed 
by means of ray-tracing simulations. With the common viscosity implementation, adiabatic gas has little effect on strong 
cluster lensing, while lower viscosity allows stronger turbulence, thus higher non-thermal pressure and a generally broader 
gas distribution, which tends to lower lensing cross sections. Conversely, cooling and star formation steepen the core density 
profiles and can thus increase the strong-lensing efficiency considerably. 



1. Introduction 

How is strong lensing by galaxy clusters affected by the intra- 
cluster gas? So far, only the dark-matter distribution was con- 
sidered for lensing, and it was ignored that on average approx- 
imately 10% to 15% of the cluster matter is contributed by the 
baryonic gas. This apparently small contribution may give rise 
to effects which can potentially change strong lensing by clus- 
ters considerably. Having seen that the central matter concen- 
tration in clusters, their asymmetries, and their substructures 
are all very important for the lensing effects they may ere- 
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in fact have substantial consequences. 

First, the finite pressure of the hot gas prevents it from 
forming density cusps. Thus, the gas distribution has to develop 
a more or less extended flat core which reduces the central den- 
sity compared to a cluster composed only of dissipation-less 
dark matter. Second, the isotropic gas pressure tends to reduce 
asymmetries in the matter distribution and to reduce the matter 
concentration in cluster substructures. This, in turn, reduces the 
gravitational tidal field (the shear) of the cluster mass distribu- 
tion. Both effects tend to act against strong cluster lensing and 
are, therefore, suspected to reduce the strong-lensing efficiency. 

On the other hand, once gas is dense enough to cool on suf- 
ficiently short time scales, it can flow towards the centre and 
accumulate there. If it forms stars instead, those act as another 
component of dissipation-less matter and may build up matter 
densities in cluster cores which exceed those expected in dark- 



matter clusters. It is not clear from the start which of these ef- 
fects will dominate and whether their competition always leads 
to the same result. As clusters form through subsequent merg- 
ers, their gas content will also modify the dynamics of the sub- 
halo accretion, shrinking the orbits of infalling sublumps and 
transporting part of their orbital angular momentum to the main 
body of the forming cluster due to friction and dissipation in 
the gas. Thus, the gas physics will change not only the density 
profile, the level of substructure, and the degree of asymmetry 
of clusters, but also modify their internal dynamics following 
merger events. 

We use gas-dynamical cluster simulations here to investi- 
gate the impact of gas physics on strong lensing by clusters. 
Most galaxy clusters were simulated five times; with adiabatic 
gas, with adiabatic gas and a novel implementation of artificial 
viscosity, with cooling gas which forms stars, with cooling gas 
which forms stars and conducts heat, and with only dark matter 
for reference. In addition, we extracted cluster-sized halos from 
a super-cluster simulation which was performed with only dark 
matter and including adiabatic gas. All simulations for a single 
cluster start from the same initial conditions at early times so as 
to allow the cluster simulations to be directly compared. This 
cluster sample is described in Sect. 2. Its strong-lensing effects 
are studied in Sect. 3, and the results are summarised and dis- 
cussed in Sect. 4. 

2. The cluster sample and its physical properties 

Our simulations were carried out with GADGET-2, a new 
version of the parallel TreeSPH simulation code GADGET 
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JSpringel etaljj 2001). It us es an en t ropy-c onserving formula- 
tion of SPH (Springel & Hernquist, 2002), and upon request 
provides radiative cooling, heating by a UV background, and a 
treatment of star formation and feedback processes. The latter 
is based on a sub-resolution model for the multi-phase stru cture 
of the interstellar medium ( Spri ngel & Hernauisl l2003ah . For 
some of the cluster simulations we also used the new method 
for describing heat conduction in SPH, which is both stable and 
manifestly conserves therma l energy even when individual and 
adaptive time-steps are used ( Jubelgas et al.ll2004T) . This imple- 
mentation assumes an isotropic effective conductivity parame- 
terised as a fixed fraction of the Spitzer rate. It also accounts 
for saturation, which can become relevant in low-density gas. 

The usual parameterisation of the a rtificial viscosity 
llMonaghan & Gmgoldlfl983llBalsaraUl995l) in SPH for an in- 
teraction of two particles a and b includes terms to account for 
shear and bulk viscosity. It is switched on only if particles are 
approaching; and for usual cosmological SPH simulations, it 
can be written as 
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Here c a t, p a b, and h a \, are the arithmetic means of the sound 
speed, the density, and the smoothing length, respectively. 
Tab - r a - Tb and v a b = V a - v a are the interparticle distance 
and the relative velocity. f a b is the mean between particles a 
and b of the viscosity limiting factor 
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which avoids spurious angular momen tum and vorticity trans- 
fer in galactic disks, as suggested by ISteinmetzl lll99d) . The 
usual choice for the parameters is a = 0.75, /3 = 2a, rj = 
0mh ab , and <r, = 0.0001c,//z,, 

Some simulations were carried out using a modified ar- 
tificial viscosity scheme suggested by iMorris &M onaghan 
Jl997h . where every particle evolves its own viscosity parame- 
ter a,, which changes with time according to 
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This causes a, to decay to a minimum value a mm = 0.01 with 
e-folding time r, which we adjust so that c,- decays over two 
smoothing lengths after the shock. The source term S ,, which 
causes a, to grow as particles approach a shock, was assumed 
to be 



S i = S*f i mzK(0,-\(V-v i ) i \~). 
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We choose S* — 0.7. Further details on this implementation 
and its consequences for the generation of turbulence within 
the intra-cluster medium will be described by Dolag et al. 2005 
(in prep). 



2.1. The cluster sample 

We used simulations of four massive galaxy clusters spanning a 
mass-range between 1.3xlO 15 /z _1 M and2.3xlO 15 /r 1 M . The 
cluster regions were extracted from a dissipation-less (dark- 
matter-only) simulation with a box-size of 479/i _l Mpc of a 
flat ACDM m odel with Q n = 0.3 , h = 0.7, cr 8 = 0.9, and 
Q b = 0.04 (see lYoshida et alfcOOlb . 

Using the "Zoomed Initial Conditions" (ZIC) technique 
(ITormenl Il997l) . they were re-simulated with higher mass and 
force resolution by populating their Lagrangian volumes in the 
initial domain with more particles, appropriately adding addi- 
tional small-scale power. The initial particle distributions (be- 
fore displacement) are of glass type dWhiteHl996l) . 

Gas was introduced into the high-resolution region by split- 
ting each parent particle into a gas and a dark-matter parti- 
cle. Thereby, the gas and the dark-matter particles were dis- 
placed by half the original mean inter-particle distance, such 
that the centre-of-mass and the momentum were conserved. 
The final mass resolution of these simulations was mdm = 
1.13 x 10 9 h~ l M e and m gas = 1.7 x 10 s hT x M e for dark mat- 
ter and gas particles within the high-resolution region, respec- 
tively. Thus, the clusters were resolved with between 2 x 10 6 
and 4x 10 6 particles, depending on their final mass. In addition, 
we used the three most massive haloes from a re-simulation 
of a super-cluster region, originating from the same cosmo- 
logical parent simulation, performed with the same resolution. 
These three haloes range in mass between 0.8 x 10 I5 /i _I M o and 
1.5 x 10 15 /z _1 M o . 

For all simulations, the gravitational softening length was 
kept fixed at e = 30.0 h kpc comoving (Plummer-equivalent), 
and was switched to a physical softening length of e = 
5.0 kpc at 1 +z = 6. 

Selection of the initial region was done with an iterative 
process involving several low-resolution, dissipation-less re- 
simulations to optimise the simulated volume. The iterative 
cleaning process ensures that all these haloes are free of con- 
taminating boundary effects up to at least 3 to 5 times the virial 
radius. With this usable volume being relatively large, the sim- 
ulation also accurately resolves the clusters' vicinity and thus 
includes the effects of all the filaments connected to the cluster. 

We used five types of simulations of this galaxy-cluster 
set. They comprise simulations with only dark matter (DM), 
simulations following the adiabatic evolution of gas but ignor- 
ing radiative cooling (GAS), and simulations including radia- 
tive cooling, heating by a UV background, and a treatment of 
the star formation and feedback processes (CSF). The feed- 
back scheme was calibrated to produce a wind velocity of 
w 350kms _I . Another simulation type we used additionally 
includes thermal conduction, for which a fixed rate of k — 1/3 
times the Spitzer rate was chosen (CSFC). This choice for k 
is appropriate in the presenc e of magnetised domains with ran- 
domly oriented fi-fields (e.g.lSaraz inil 19881 ) or for a ch aotically 
tangled magnetic field (Naravan & Medvedev, 2001). Finally, 
we used one type of adiabatic gas simulations in which we ap- 
plied an implementation of the artificial viscosity, which damps 
the build-up of viscosity in the time domain and thus reduces 
it considerably where it is not needed (GAS_NV). The lower 
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artificial viscosity and the absence of a limiting physical vis- 
cosity allow strong turbulence to build up in the centres of the 
simulated clusters. Its contribution to the pressure leads to a 
significant density decrease in the cluster cores. For more de- 
tails, see Dolag et al. 2005 (in prep.). 



2.2. Halo shape and particle angular momentum 

Before discussing the lensing properties of our simulated 
galaxy clusters, we compare their shapes, density profiles, and 
angular momentum distributions for the five different physical 
gas models used in the simulations. Figure^shows typical pro- 
files of the total density (dark matter and baryons) for the dif- 
ferent types of gas physics. The most obvious difference is the 
steeper inner slope in the simulations with cooling and star for- 
mation. Although we use a state-of-the -art impl ementati on of 
cooling , feedback, and star formation ( Spring e! & HernquisA 
l2003bl) we should point out that it is not entirely clear how 
realistic the profiles for these simulations are close to the clus- 
ter centre, as the central cD galaxy will contribu te substan tially 
to the core density profiles (Le wis et all 12000: Yos hida et all 
I2OO2I) . The core density of stars in the simulations is larger than 
observed (the simulated clusters seem to suffer t o some degree 
from over-cooling), but some authors (see lKravtsov et alJ2 005 ) 
argue that part of the discrepancy may be due to stellar mass 
missed in the observations. 

Despite the isotropic thermal gas pressure, the density pro- 
file in the GAS model is not significantly shallower than in the 
DM simulation. This can be understood from the fact that gas 
particles can reduce their angular momentum by collisions (see 
below), which lets them sink towards the cluster centre more 
easily. In the GAS_NV simulation, the additional pressure sup- 
port due to strong turbulence allowed by the lower viscosity re- 
duces the gas density in the inner region of the simulated cluster 
(see Dolag et al. 2005, in prep., for more detail). However, the 
impact of the turbulence on the density profiles of real clusters 
requires further investigation because there the physical viscos- 
ity, which is not yet included in the simulations, may or may not 
limit the amount of turbulence to smaller values. 

The effects discussed here can also be seen in Fig. [2] which 
shows the baryonic mass fractions for the GAS and GAS_NV 
simulations, and the fractions of gas, stars, and the total baryon 
fraction for the CSF simulation of cluster gl as a function of 
the radius of the sphere around the cluster centre in which it 
was computed. The baryon fractions of the GAS and CSF sim- 
ulati ons are also in goo d agreement with the results obtained 
bv lKravtsov et alJ J2005I) . 

Figure [5] displays the mass fraction of the particles with 
specific angular momentum \L\/m < j against the threshold 
j for the cluster gl at redshift z = 0.2975. Here, L and m 
are the angular momentum and the mass of the cluster parti- 
cles, respectively. This quantity is plotted for the particles of 
the simulations with only dark matter and for the gas and dark 
matter particles of the simulation including adiabatic gas. The 
angular-momentum profiles of the GAS_NV, CSF, and CSFC 
simulations, which are not plotted for clarity, are qualitatively 
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Fig. 1. Profiles of the total density of cluster gl at redshift z 
0.2975 for different gas-physical models. 
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Fig. 2. Baryonic mass fraction M <r /M tota i <r in the GAS and 
GAS_NV simulations, and mass fractions of gas, stars, and all 
baryons in the CSF simulation, of cluster gl at redshift z = 
0.2975 as a function of distance r from the cluster centre. 



similar to the GAS case. Only particles with a distance smaller 
than 250 h kpc from the cluster centre were included. 

We see from this figure that the profiles of the dark-matter 
particles in the two different simulations are almost identical. 
However, the gas particles typically have a significantly lower 
specific angular momentum, and the same behaviour is also 
found in the simulations with cooling and star formation and 
in the simulations with the new model for artificial viscosity. 
When studying the time evolution of these profiles, we found 
that the specific angular momenta of both dark matter and gas 
particles are boosted towards higher values during mergers. 
Afterwards, the angular-momentum profile of the dark mat- 
ter is almost conserved, while the gas relaxes and the parti- 
cles lose specific angular momentum in collisions. This hap- 
pens because the cluster halos lack a well-defined rotation axis 
and the orbital planes of gas particles have essentially random 
orientations. Collisions tend to average out differences in orbit 
orientation and thereby reduce the specific angular momentum 
of gas particles. Therefore, the difference between the specific 
angular momentum of dark matter and gas increases with each 
merger. We thus find somewhat larger deviations in large halos, 
which have on average experienced more mergers. Reducing 
their specific angular momenta allows the gas particles to sink 
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more easily towards the cluster centre. A similar behaviour 
of the specifi c angular momenta of gas and dark matter was 
found bv lNavarro & Ben j dl99lh and iNavarro & White! (1 1994 
in galaxy-sized halos. 

By comparing the monopole to the quadrupole moment in 
circular she lls around the cluster cen tre, we recover the result 
obtained bv lKazantzidis et al](l2004l) . who found that halos are 
more spherical in the simulations with cooling, feedback, and 
star formation than in the dissipation-less and adiabatic gas 
simulations. For example at the scale radius we find 15 percent 
smaller ratios of quadrupole to monopole moments in projec- 
tions of simulations with cooling and star formation. We also 
find more substructure in the form of small clumps in the CSF 
and CSFC models. 
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Fig. 3. Cumulative distributions of the specific angular momen- 
tum of the particles contained in the central region of cluster gl 
at z — 0.2975, within a radius of r — 250 h~ x kpc. 



3. Strong lensing cross sections 

3.1. Numerical methods 



From the projected mass map, we calculated the conver- 
gence k and its Fourier transform k, which is related to the 
Fourier transform of the lensing potential {j/ by 



(6) 



Employing fast-Fourier methods for deriving the lensing po- 
tential according to this equation is, however, problematic. 
Discrete Fourier transform algorithms assume that the func- 
tion to be transformed is periodic on its support, which is 
not the case for isolated and finite cluster convergence fields. 
Equation 0, therefore, does not yield the lensing potential of 
a single cluster, but that of an infinite, two-dimensional array 
of clusters (see Fig.@J, in which the original cluster is repeated 
on a grid whose periodicity is set by the side length of a single 
cluster field (6 h' 1 Mpc comoving in our case). 

Thus, a sufficiently accurate result is only achieved close 
to the centre of the cluster field. For reducing the error, we 
could surround the cluster by an even larger zero-padded field. 
However, doing so without losing resolution substantially in- 
creases the demands on computer memory and slows down the 
computation. We have therefore chosen to use a new method to 
correct for these errors. We place a point mass that concentrates 
the total projected mass of the cluster at the cluster's centre-of- 
mass. We then calculate the lensing potential for an array of 
these point masses with the Fourier method mentioned above. 
Next, we analytically subtract the lensing potential of the single 
point mass located at the cluster's centre-of-mass so as to ob- 
tain the lensing potential of an array of point masses with one 
mass missing at the position of the original cluster. Subtracting 
this from the potential obtained for the array of clusters, we 
correct for the additional clusters implicitly produced by the 
fast-Fourier algorithms, which assume a periodic convergence 
array. The remaining error comes only from the higher multi- 
pole moments of these additional clusters and can be neglected 
in the central sixteenth (1024 x 1024 points) of the grid which 
we use for doing arc statistics. From this corrected lensing po- 
tential, we calculate the deflection angle. 



The centre of each clu ster is found by usi ng the halo-finder 
algorithm discussed in Torm en et al . (2004). It estimates the 
dark-matter density at the position of each dark-matter particle 
by determining the distance to the tenth-closest neighbour dio 
and by also assuming that the density at the particle position is 
proportional to djj? . Starting there, the virial sphere of the parti- 
cle distribution is found, in which the gravitational potential is 
determined. The halo centre is then taken to be at the potential 
minimum. 

For studying the lensing properties of each cluster, we 
chose a sphere of comoving radius 3 hr x Mpc around the cluster 
centre and projected all cluster particles inside this region onto 
an equidistant grid with a resolution of 4096 x 4096 cells. For 
the projection we used the same cubic spline function which 
is used in the GADG ET code as the SPH smoothing kernel 
(Snrin gel et all hoOl). We projected by calculating the overlap 
of the projected spline function and the square representing the 
pixels of the grid. 




cluster 

field for calculation of lensing 
potential 

"additional cluster" from assuming 
periodicity 

point mass to correct for 
"additional cluster" 
field used for arc statistics 



Fig. 4. Illustration of our method for correcting the error caused 
by the implicit assumption of a periodic input function in fast- 
Fourier techniques. 



For finding the images of a number of sources large enough 
for statistical analysis, w e follow the method introduced by 
iMiralda-Escudd (TT993a b) and adapted to non-analytic models 
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bv iBartelmann & Weissl ( 1 1994 and iBartelmann et all dl995h . 
A previo us version of our algorithm is also discussed in some 
detail in iMeneehetti etail §000). It places elliptical sources 
with an equivalent radius of 0.5 arc seconds on an adaptive 
grid in the source plane, which is fixed at redshift 1.5 (the lens 
redshift is taken to be the redshift of the simulation snapshot), 
such that there is an increasing density of sources close to the 
caustics. A statistical weight is assigned to each source, which 
is given by the area represented by the source. 

Next, the deflection angles are used to trace light rays back- 
wards and map each grid point from the lens plane to the source 
plane. The images of a source are found by checking which grid 
points, when mapped back to the source plane, are enclosed by 
the ellipse corresponding to the source considered. To deter- 
mine the image properties (e.g. length L, width W, and curva- 
ture radius R), we find the image point (a) which, when mapped 
to the source plane, falls closest to the source centre. We then 
find the image point (b) which is the farthest from (a), and the 
image point (c), which is the farthest from (b). The method is 
illustrated in Fig.|5] 

We could fit a circle through these three points and use the 
arc length from (b) to (c) as the length of the image. We would 
then determine the image perimeter by walking along the or- 
dered boundary points and summing up their mutual distances. 
But since grid cells in the lens plane are only classified as be- 
longing to the image if their centres fall within the source, the 
boundary points of the image (including (b) and (c)) will on 
average be about half a grid constant further inside the image 
than the true perimeter. Thus, we would systematically under- 
estimate the length of the image by roughly one grid constant. 
We would also underestimate the perimeter of the image. We 
correct for this by adding one grid constant to the arc length 
from (b) to (c) to obtain the image length, and four grid con- 
stants to the sum of the distances of the boundary points to find 
the image perimeter. This is also shown in Fig. [5] We found 
that these new corrections further reduce the weak dependence 
of the lensing cross section on the resolution of the grid used 
for its computation. 

The image area is calculated directly from the number of 
image points. As discussed in Meneghett i et al . (2000), we 
search a simple geometric figure (ellipse, circle, rectangle or 
ring) with equal area and length to determine the image width 
W, which is approximated by the minor axis of the ellipse, the 
diameter of the circle, the smaller side of the rectangle or the 
width of the ring, respectively. We choose the type of the figure 
by comparing its circumference to the perimeter of the image, 
which we find using the method discussed above. 

We finally determine the lensing cross section C7.5 by sum- 
ming up the statistical weights of the sources having images 
with a length-to-width ratio L/W > 7.5 and calculating the co- 
moving area corresponding to it. If there is more than one such 
image for a source, we multiply the statistical weight of this 
source by the number of these images. 




Fig. 5. Correction of image length and perimeter. 
3.2. Results 

We used the method discussed above to find the strong-lensing 
cross sections of the clusters gl, g8, g51, g72, and of the three 
largest halos in the super-cluster simulation g696. We did this 
for 43 simulation snapshots with redshifts ranging between 0. 1 
and 1.05, and for the five different physical gas models dis- 
cussed in Sect. |2 except for g696, for which only DM and 
GAS simulations are available. For each halo, we used three 
different projections, namely along the x, y and z axes of the 
simulation boxes in which the clusters are randomly oriented. 

In Fig.|S] we plot the cross section as a function of redshift 
for one of the projections of g51. We can see that despite the 
gas pressure, adiabatic gas with the standard artifical viscos- 
ity does not reduce the strong lensing cross section compared 
to the dark-matter-only simulation. Using the new artifical vis- 
cosity scheme, however, leads to somewhat smaller cross sec- 
tions due to the more extended gas distribution it implies. On 
the other hand, cooling, star formation and feedback make the 
simulated cluster a significantly more efficient lens. These three 
properties are typical for most of the clusters we have stud- 
ied. There is generally no large difference between the DM and 
(adiabatic) GAS models, a somewhat smaller cross section in 
the GAS_NV model, and cross sections larger by a factor of 
1.5 to 3 in the CSF and CSFC models. In some cases, however, 
even the adiabatic gas with standard artifical viscosity causes 
an increase in the cross section compared to the simulations 
containing only dark matter. This is illustrated in Fig.0 

For determining the impact of cluster ellipticity and sub- 
structure on the strong lensing cross sections, we transform the 
maps of the surface mass density of our cluster halos to polar 
coordinates (centred on the cluster halo) and average over the 
polar angle to obtain maps of an axially symmetrised cluster. 
Then, we use the same methods as before to compute strong- 
lensing cross sections. 

Figure [8] shows the cross sections of the GAS and CSF 
versions of g51 and of its axially-symmetrised GAS and CSF 
counterparts. One can clearly see that substructure and elliptic- 
ity significantl y increase the strong-lensin g cross section s (see 
also IBartelmann et aDll995l and Meneghet ti et alJl2003l) . Note 
also that the increase of the cross section in the CSF simula- 
tion, compared to the GAS simulation, is almost the same for 
the original cluster and its axially-symmetric variant. This can 
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Fig. 6. Strong lensing cross sections for the projection along the 
y axis of the DM, GAS, GASJSfV, CSF, and CSFC simulations 
of cluster g51, for an arc length-to- width ratio of 7.5 or more. 
The units are comoving (Mpc//j) 2 . 
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GAS simulation 
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CSF simulation 
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Fig. 7. Strong lensing cross section for the DM, GAS, 
GAS_NV, CSF, and CSFC simulation of cluster gl and an arc 
length-to-width ratio of > 7.5. 



be interpreted such that this increase in the cross section in the 
cooling and star formation simulations is caused mainly by the 
steeper density profile and not by any changes of the ellipticity 
or the substructure. 

There is also no qualitative difference between the lensing 
properties of the simulated clusters and the largest halos of the 
simulated super cluster region. 

The different physics in the simulations with only dark mat- 
ter particles compared to simulations including gas directs sub- 
halos passing close to the main halo during a merger into dif- 
ferent orbits. Compared to the dissipation-less simulation with 
only dark matter, the sub-halo loses angular momentum and 
energy in the gas-dynamical simulations and is directed into 
a less elliptical orbit. It therefore returns earlier for the next 
passage of the main halo's centre. This is illustrated in Fig. [5] 
which shows the position of a sub-halo in the rest frame of 
the main halo of the cluster g72. Note that the x and y axes are 
scaled differently for clarity. The sub-halos in the DM and GAS 
simulations initially move approximately synchronously and 
approach the main halo with almost the same velocity. Later, 
however, the different dynamics makes the GAS sub-halo's or- 



GAS sim. 
GAS sim. symmetrized 

CSF sim. 
CSF sim. symmetrized 




Fig. 8. Strong lensing cross sections for the GAS, CSF, and 
axialiy-symmetrised GAS and CSF versions of cluster g51. 



bit substantially less elliptical, hence it returns earlier for the 
next core passage. 

This effect also has an impact on the strong lensing cross 
section, because the cross section of a halo increases when 
a sub-halo cros ses the cluster core during a merger (see 
iTorri etalJEooi and iMeneghetti et alJ 12004 . In Fig. [K)| the 
peaks in the cross section corresponding to the three successive 
passages, illustrated in Fig. [9] are marked by arrows for both the 
DM and the GAS simulations. The first peak happens in both 
cases at a redshift z ~ 0.62. The second and the third peaks, 
however, occur at significantly earlier times in the GAS sim- 
ulation. The peak positions for the GAS_NV, CSF, and CSFC 
simulations are typically very similar to the GAS case. Note 
that the different heights of the first peak in the different simu- 
lations do not imply a fundamental difference, because the am- 
plitude of sharp peaks depends strongly on the exact time when 
the snapshot was taken (the time resolution is just a few snap- 
shots for the passage). Thus, a slight deviation in the timing 
may result in peaks with apparently very different heights. 

4. Summary and discussion 

We have compared the density profiles, shapes, and strong 
lensing properties of numerically simulated galaxy clusters. 
The simulations were performed starting from the same initial 
conditions, but with different gas physics. Five different gas- 
physical models were employed. They contain: 

- DM model: dark matter only; 

- GAS model: dark matter and adiabatic gas; 

- GAS_NV model: dark matter, adiabatic gas, and a new im- 
plementation of the artificial viscosity, reducing the viscos- 
ity where it is not needed; 

- CSF model: dark matter, cooling gas, a star formation 
model, and feedback; 

- CSFC model: dark matter, cooling gas, a star formation 
model, feedback, and thermal conductivity; 

Our cluster sample consisted of four simulated galaxy clus- 
ters and the three largest halos of a simulation of a super-cluster 
region. We used 43 snapshots of these halos between redshifts 
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Fig. 9. Position (in the main halo's rest frame and in comoving 
coordinates) and radial disctance from the main halo of the sub- 
halo whose passage produces the peaks in cluster g72's lensing 
cross section for the DM and GAS case. 
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0.4 0.5 
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Fig. 10. Strong-lensing cross section for the DM and GAS sim- 
ulations of cluster g72 and a length-to-width ratio of 7.5 or 
more. The arrows mark the three successive passages of the 
same sub-halo for both models. 

1 .05 and 0.1. For studying strong lensing, we used the three dif- 
ferent projections along the coordinate axes of the simulation 
volume, and sources fixed at redshift 1.5. 

We find significantly steeper inner slopes for the density 
profiles of halos simulated with cooling and star formation. On 
the other hand, adiabatic gas with a standard artifical viscosity, 
in spite of its isotropic thermal pressure, does not make density 
profiles shallower compared to the dissipation-less dark-matter 
simulations. This can be understood from the fact that gas can 
reduce its specific angular momentum by collisions. The or- 
bits of gas particles in a cluster halo are randomly oriented. 
Collisions tend to average out these differences and thereby re- 
duce the specific angular momentum of the gas, which helps 
the gas to move towards the cluster centre. This effect can also 
be seen in the specific angular-momentum profiles we present 
in Fig. [3] However, the additional pressure caused by strong 



turbulence in the GAS_NV simulation somewhat reduces the 
density close to the cluster centre. 

We then performed ray-tracing simulations for the numer- 
ically simulated galaxy clusters to study their strong lensing 
properties. We calculated cross sections for long thin arcs with 
a length-to- width ratio equal to or larger than 7.5. For the simu- 
lations with cooling and star formation, we found significantly 
larger cross sections. Thermal conductivity has no big impact 
on strong lensing. Despite its pressure, adiabatic gas with stan- 
dard artifical viscosity does not reduce the cross section for 
long thin arcs. In some cases, it even makes the cluster a more 
efficient lens. On the other hand, simulating adiabatic gas with 
the new scheme for reduced artifical viscosity leads to some- 
what smaller strong lensing cross sections compared to the DM 
and GAS runs. 

Note that despite our use of a state-of-the-art model for 
cooling, star formation, and feedback, the simulated clusters 
suffer from some over-cooling. Therefore, the density close to 
the centre and the increase of the strong lensing cross section 
may be over-predicted in the simulations with cooling and star 
formation compared to real clusters. Also, the impact of tur- 
bulence on galaxy clusters needs further investigation. In real 
clusters, the physical viscosity of the cluster gas, which is not 
yet included in simulations, will regulate turbulence and may 
lead to a smaller amount of turbulence than in the GAS_NV 
simulations. It will also be interesting to investigate the effect 
of this new scheme of artifical viscosity on simulations includ- 
ing cooling, star formation, and feedback. 

To study the impact of ellipticity and substructure on the 
lensing properties we transformed maps of the surface mass 
density of our cluster halos to polar coordinates and averaged 
over the angle, so that we obtained maps of an azimuthally 
symmetrised halo. We compared the cross sections of these 
axially symmetrised clusters to those of the numerically sim- 
ulated ones and found that ellipticity and substructure are im- 
portant for strong lensing in all of the five gas physical models. 
However, their impact on the cross section turns out to be very 
similar in the simulations without and with cooling and star 
formation. Thus, the larger cross sections we obtained for sim- 
ulations with cooling and star formation are mainly caused by 
the steeper density profile and not by changes in ellipticity and 
substructure. 

Sub-halos passing close to the main halo lose angular mo- 
mentum and energy and are directed into less elliptical orbits 
in simulations with gas compared to dark-matter-only simula- 
tions. During a merger, such a sub-halo will return at an earlier 
time for the next passage. Thus, mergers proceed faster in sim- 
ulations with gas. This can also be seen in the strong lensing 
cross sections: typically, the peak corresponding to the first pas- 
sage of a sub-halo happens at the same time in the dissipation- 
less and the gas-dynamical simulations, but the peaks corre- 
sponding to the next passages are shifted to earlier times in the 
simulations with gas. 

The work presented here clearly shows that cluster gas 
can have a significant impact on strong lensing cross sections. 
However, the importance of the different effects (turbulence, 
cooling, star formation, mergers) is different for each individ- 
ual cluster and changes during the cluster's evolution. To infer 
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cluster properties from observations and for studying the im- White, S. 1996, in Cosmology and Large-Scale Structure, ed. 
portance of these effects in real clusters, it is thus more promis- R- Schaeffer, J. Silk, M. Spiro, & J. Zinn- Justin (Elsevier, 
ing to model observed clusters individually than to do statistical Dordrecht), 395 

analyses of large cluster samples. For this it will be interest- Yoshida, N, Sheth, R., & Diaferio, A. 2001, MNRAS, 328, 669 
ing to combine the lensing information with thermal Sunyaev- Yoshida ' N - Stoehr ' F > S P rin g e1 ' v - & White ' S ' 2002 ' MNRA S, 335, 
Zeldovich effect and X-ray observations. We are planning to 
simulate X-ray emission and thermal SZ effect for the numeri- 
cal galaxy clusters used here and to study methods for combin- 
ing such (mock) observations to model individual clusters and 
reconstruct their mass distributions. 
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